%# Field 1: fday
%# Field 2: longitude (degrees)
%# Field 3: latitude (degrees)
%# Field 4: local time (hours)
%# Field 5: season (day of year)
%# Field 6: eastward electric field (V/m)
%# Field 7: DC current shift (A/m)
%# Field 8: CHAMP peak current value (A/m)
%# Field 9: KP
%# Field 10: F10.7 (W/m^2)
%# Field 11: F10.7A (W/m^2)
%# Field 12: R^2 (coefficient of determination)
%# Field 13: R (correlation of CHAMP and model profiles)
%# Field 14: chi^2


for i = 1:length(eef),
    if eef(i,4) > 7.0 & eef(i,4) < 17.0,
fprintf(fp,'%10.6f %8.3f %10.6f %10.6f %10.6f\n',eef(i,1),eef(i,2),eef(i,4),eef(i,5),eef(i,11));
    end;
end;